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Abstract 

We show that it is possible to improve the infrared aspects of the standard treatment of 
the DGLAP evolution theory to take into account a large class of higher order corrections 
that significantly improve the precision of the theory for any given level of fixed-order 
calculation of its respective kernels. We illustrate the size of the effects we resum using 
the moments of the parton distributions. 
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In the preparation of the physics for the precision QCDxEW(electroweak) [1,2] LHC 
physics studies, all aspects of the calculation of the cross sections and distributions for 
the would-be physical observables must be re-examined if precision tags such as that 
envisioned for the luminosity theoretical precision are to be realized, i.e., 1% cross section 
predictions for single heavy gauge boson production in 14 TeV pp collisions when that 
heavy gauge boson decays into a light lepton pair. The QCD DGLAP [3] evolution of 
the structure functions from the typical reference scale of data input, /io ~ 1 — 2GeV, 
to the respective hard scale is one step that warrants further study, as it is well-known 
to many. Many authors [4-7] have provided excellent realizations of this evolution in 
the recent literature. Here, we will re-examine the infrared aspects of the basic DGLAP 
theory itself to try to improve the treatment to a level consistent with the new era of 
precision QCDxEW physics needed for the LHC physics objectives. 

Specifically, the motivation for the improvement which we develop can be seen already 
in the basic results in Refs. [3] for the kernels that determine the evolution of the structure 
functions by the attendant DGLAP evolution of the corresponding parton densities by 
the standard methodology. Consider the evolution of the non-singlet (NS) parton density 
function q NS (x), where x can be identified as Bjorken's variable as usual. The basic 
starting point of our analysis is the infrared divergence in the kernel that determines this 
evolution: 

dq NS (x,t) a s (t) fUy 



dt 2?r 



/ % NS (y,t)P qq (x/y) (1) 



where the well-known result for the kernel P qq (z) is, for z < 1, 

PM = C F 1 -^- (2) 

when we set t = ln/i 2 //^ for some reference scale /io with which we study evolution to 
the scale of interest /i. Q| Here, Cp = (iV 2 — 1)/(2N C ) is the quark color representation's 
quadratic Casimir invariant where N c is the number of colors and so that it is just 3. This 
kernel has an unintegrable IR singularity at z — 1, which is the point of zero energy gluon 
emission and this is as it should be. The standard treatment of this very physical effect 
is to regularize it by the replacement 

(3) 



(1-2) (1-*) + 

with the distribution ^ z ^ + defined so that for any suitable test function f(z) we have 

l dz M = /yw-/w. (4) 







z )+ Jo (! - z ) 



1 We will generally follow Ref. [8] and set /io = h-QCD without loss of content since dt = dt' when 
t = In /i 2 / Aq CD , t' = hi/i 2 //ig for fixed values of Aqcd,^o- 
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A possible representation of 1/(1 — z)+ is seen to be 



1-z 



l-e-z)+lne6(l-z) (5) 



with the understanding that e j 0. We use the notation 8(x) for the step function from 
for x < to 1 for x > and 5(x) is Dirac's delta function. The final result for P qq {z) is 
then obtained by imposing the physical requirement [3] that 

dzP qq (z) = 0, (6) 



o 



which is satisfied by adding the effects of virtual corrections at z = 1 so that finally 

P,,(z) = C F (l±^y-6 ( l- Z )). (7) 



The smooth behavior in the original real emission result from the Feynman rules, 
with a divergent 1/(1 — z) behavior as z —* 1, has been replaced with a mathematical 
artifact: the regime 1 — e < z < 1 now has no probability at all and at z = 1 we have a 
large negative integrable contribution so that we end-up finally with a finite (zero) value 
for the total integral of P qq (z). This mathematical artifact is what we wish to improve 
here; for, in the precision studies of Z physics [9-11] at LEP1, it has been found that such 
mathematical artifacts can indeed impair the precision tag which one can achieve with 
a given fixed order of perturbation theory. An analogous case is now well-known in the 
theory of QCD higher order corrections, where the FNAL data on spectra clearly show 
the need for improvement of fixed-order results by resumming large logs associated with 
soft gluons [12,13]. For reference, note that at the LHC, 2 TeV partons are realistic so 
that z = 0.001 means ~ 2 — 3 GeV soft gluons, which are clearly above the LHC detector 
thresholds, in complete analogy with the situation at LEP where z = 0.001 meant ~ 100 
MeV photons which were also above the LEP detector thresholds - just as resummation 
was necessary to describe this view of the LEP data, so too we may argue it will be 
necessary to describe the LHC data on the corresponding view. And, more importantly, 
why should we have to set P qq {z) to for 1 — e < z < 1 when it actually has its largest 
values in this very regime? 

By mathematical artifact we do not mean that there is an error in the computations 
that lead to it; indeed, it is well-known that this +- function behavior is exactly what one 
gets at 0(a s ) for the bremsstrahlung process. The artifact is that the behavior of the 
differential spectrum of the process for z — > 1 in 0(a s ) is unintegrable and has to be 
cut-off and thus this spectrum is only poorly represented by the 0(a s ) calculation; for, 
the resummation of the large soft higher order effects as we present below changes the 
z — > 1 behavior non-trivially, as from our resummation we will find that the ^-^-behavior 
is modified to (1 — z)' y ~ 1 ,'-f > 0. This is a testable effect, as we have seen in its QED 
analogs in Z physics at LEP1 [9-11]: if the experimentalist measures the cross section 
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for bremsstrahlung for gluons (photons) down to energy fraction eg, €q > 0, in our new 
resummed theory presented below, the result will approach a finite value from below as 
— eJJ whereas the 0(a s ) +-function prediction would increase without limit as — lneo- The 
exponentiated result has been verified by the data at LEP1. 

The important point is that the traditional resummations in N-moment space for 
the DGLAP kernels address only the short- distance contributions to their higher order 
corrections. The deep question we deal with in this paper concerns, then, how much of the 
complete soft limit of the DGLAP kernels is contained in the anomalous dimensions of the 
leading twist operators in Wilson's expansion, an expansion which resides on the very tip 
of the light-cone? Are all of the effects of the very soft gluon emission, involving, as they 
most certainly do, arbitrarily long wavelength quanta, representable by the physics at the 
tip of the light-cone? The Heisenberg uncertainty principle surely tells us that answer 
can not be affirmative. In this paper, we calculate these long-wavelength gluon effects 
on the DGLAP kernels that are not included (see the discussion below) in the standard 
treatment of Wilson's expansion. We therefore do not contradict the results of the large 
N-moment space resummations such as that presented in Ref. [14] nor do we contradict 
the renormalon chain-type resummation as done in Ref. [15]. 

We employ the exact re-arrangement of the Feynman series for QCD as it has been 
shown in Ref. [16, 17]. For completeness, as this QCD exponentiation theory is not gen- 
erally familiar, we reproduce its essential aspects in our Appendix. The idea is to sum 
up the leading IR terms in the corrections to P qq with the goal that they will render 
integrable the IR singularity that we have in its lowest order form. This will remove the 
need for mathematical artifacts and exhibit more accurately the true predictions of the 
full QCD theory in terms of fully physical results. 



We apply the QCD exponentiation master formula in eq. flHOj) in our Appendix (see 
also Ref. [16]), following the analogous discussion then for QED in Refs. [10,11], to the 
gluon emission transition that corresponds to P qq (z), i.e., to the squared amplitude for 
q — > q(z) + G(l — z) so that in the Appendix one replaces everywhere the squared ampli- 
tudes for the Q'Q — > Q"'Q" processes with those for the former one plus its nG analoga 
with the attendant changes in the phase space and kinematics dictated by the standard 
methods; this implies that in eq.(53) of the first paper in Ref. [3] we have from eq. ffSOl) 
the replacement ( see Fig. ) 



p BA = p° A = Ui-')Y. 1 
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2 ±r. n 



spins 
\V A 



Pba = L ~z{l ~z)J2 ^ A ^ 2 z^F YFS { lq )e^ 
£ p I 



spins 



where A = q, B = G, C = q and Va^b+c is the lowest order amplitude for q 
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(a) (b) 
Figure 1: In (a), we show the usual process q — > q(l — z) + G(z); in (b), we show its 
multiple gluon improvement q — > q(l — z) + Gi(£i) + • ■ ■ + G n (^ n ), z = ^ . ^. 



+ q(l — z), so that we get the un- normalized exponentiated result 



where [10,11,16,17] 



P qq {z) = C F F YFS { lq )e^ l -^{l - zY* (9) 



a s AC F 

7 » = Cf ^ = a" (10) 

4, = a + _ J, (ii) 



2 vr v 3 2 

and 



W7.) = f^y (12) 

Here, 

2 

A) = 11 - gW/ 

, where n/ is the number of active quark flavors, 

C E = .5772 . . . 

is Euler's constant and T(w) is Euler's gamma function. The function Fyfs{ z ) was already 
introduced by Yennie, Frautschi and Suura [18] in their analysis of the IR behavior of QED. 
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We see immediately that the exponentiation has removed the unintegrable IR divergence 
at z — 1. For reference, we note that we have in ([9]) resummed the terms^l 0(ln fc (l — 
z)t e a™), n > £ > k, which originate in the IR regime and which exponentiate. The 
important point is that we have not dropped outright the terms that do not exponentiate 
but have organized them into the residuals j3 m in the analog of eq.( l80|) . The application 
of eq. ([80T) to obtain eq.® proceeds as follows. First, the exponent in the exponential 
factor in front of the expression on the RHS of eq. (180]) is readily seen to be from eq. (!77|) . 
using the well-known results for the respective real and virtual infrared functions from 
Refs. [16,17], 

SUM IR (QCD) = 2a s ReB QCD + 2a s B QCD {K max ) 

1 fnn Kmax , n as + ■ a sC F/ T^ 2 1 N \ ( 13 ) 

— | ZUp — tin — — \- Lsf — t 



2 \ 7r E 2ir ix v 3 2 

where on the RHS of the last result we have already applied the DGLAP synthesization 
procedure in the third paper in Ref. [16] to remove the collinear singularities, In Aq CD / mi — 
1, in accordance with the standard QCD factorization theorems [19]. This means that, 
identifying the LHS of eq. (1801) as the sum over final states and average over initial states 
of the respective process divided by the incident flux and replacing that incident flux 
by the respective initial state density according to the standard methods for the pro- 
cess q — > q(l — z) + G(z), occurring in the context of a hard scattering at scale Q as 
it is for eq.(53) in the first paper in Ref. [3], the soft gluon effects for energy fraction 

< z = K max /E give the result, from eq. flHOl) . that, working through to the level and 
using q2 to represent the momentum conservation via the other degrees of freedom for the 
attendant hard process, 

J p BAdtd Z = e SUM IR (QCD) (z ) J0 q J ^^fe^H/^^^W^'-l]} 

rf3fc lS/;„N f ^iyipi-P2-k 1 )+f k<Kmax ^S QC D(k)[e-^- k -l]} 



d 3 p 2 d 3 q 2 

' ' ' ' J o n 



P2 <?2 

SUM IR (QCD)(z) [r% f°° d V fi y .(E 1 -E 2 )+f k<Kma * 4^ QC d W[e-'^-l]} 

(2tt) 



e SUM IR ^UUj(zj J {/3 
3k. _ poo 



oo 



d ^Mki) r — c {iyiEl - E2 - k ° )+ ^ <KmaX ^gQCD(fc)[e- itf - fc -l]> 



1 d3 P2 

% 2 V 



(14) 
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Following the standard LEP Yellow Book [9] convention, we do not include the order of the first 
nonzero term in counting the order of its higher order corrections. 
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where we set Ei — pf, i — 1,2 and the real infrared function Sqcd(&) is well-known as 
well: 

5 n\ _ a » C F /Pi P2V , 1K x 

^QCDlfcJ I ^ - ^ I IDGLAP synthesized 

and we indicate as above that the DGLAP synthesization procedure in Refs. [16] is to be 
applied to its evaluation to remove its collinear singularities; we are using the kinematics 
of the first paper in Ref. [3] in their computation of Pba(z) in their eq.(53), so that 
the relevant value of k\ is indeed Q 2 . It means that the computation can also be seen 
to correspond to computing the IR function for the standard t-channel kinematics and 
taking | of the result to match the single line emission in Pc q - The two important integrals 
needed in ( IT4|) were already studied in Ref. [18]: 



I YFS (zE,0) = I ^-e 



DC 



I Y fs(zE, fci) = / — e 



(16) 



7F, 

f-^IypsizE^) 



zE - ki 

When we introduce the results in ([TBI) into (I14p we can identify the factor 

/ + / ^M^r^) || = / *S^PL<b + 0( oS. 

(17) 

where Pg^ is the unexponentiated result in the first line of ([8]). This leads us finally to 
the exponentiated result in the second line of (jSj) by elementary differentiation: 

P BA = P BA z^F YF s{l q )e^" (18) 

Here, we also may note how one can see that the terms we exponentiate are not 
included in the standard treatment of Wilson's expansion: From the standard method [20], 
the N-th moment of the invariants T^g, i — L, 2, 3, I = q, G, of the forward Compton 
amplitude in DIS is projected by 



V 



N 



Nl dp^ - --dp 



\p=0 



(19) 



where x B j = Q 2 /{2qp) in the standard DIS notation; this projects the coefficient of 
1/(23^)^. For the dominant terms which we resum here, the characteristic behavior 
would correspond formally to 7 g -dependent anomalous dimensions associated with the 
respective coefficient whereas by definition Wilson's expansion does not contain such. In 
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more phenomenologically familiar language, it is well-known that the parton model used 
in this paper to calculate the large distance effects that improve the kernels contains 
such effects whereas Wilson's expansion does not: for example, the parton model can 
be used for Drell-Yan processes, Wilson's expansion can not. Simalarly, any Wilson- 
expansion guided procedure used to infer the kernels via inverse Mellin transformation, 
by calculating the coefficient of (l/z) n in Wilson's expansion, will necessarily omit the 
dominant IR terms which we resum. Here, we stress that we refer to the properties of the 
expansion of the invariant functions Tj, not to the expansion of the kernels themselves, as 
the latter are related to the respective anomalous dimension matrix elements by inverse 
Mellin transformations. 



The normalization condition in eq.([6]) then gives us the final expression 

'l + z 2 



P qq (z) = C F F YFS { lq )e* 6 « 



I - z 



(f - z)** - fM5(l - z) 



(20) 



where 



'■<*> = k " ^TT + ^2- (21) 

The latter result is then our IR-improved kernel for NS DGLAP evolution in QCD. We 
note that the appearance of the integrable function (1 — z)~ l+lq in the place of ^j 1 ^ 
was already anticipated by Gribov and Lipatov in Refs. [3]. Here, we have calculated the 
value of 7 9 in a systematic rearrangement of the QCD perturbation theory that allows 
one to work to any exact order in the theory without dropping any part of the theory's 
perturbation series. 

The standard DGLAP theory tells us that the kernel Pcq^z) is related to P qq (l — z) 
directly: for z < 1, we have 

P Gq (z) = P qq (l -z)= C F F YFS { lq )e^ 1 + ^~ z) z*. (22) 

This then brings us to our first non-trivial check of the new IR-improved theory; for, the 
conservation of momentum tells us that 



/ dzz(P Gq (z)+P qq (z)) = 0. 
Jo 



(23) 



Using the new results in eqs. ( !20|22l) . we have to check that the following integral vanishes: 
/ = f dzz + + i±^(l - zf« - f q { lq )5{l - z)) . (24) 



o 



z 

To see that it does, note that 



j-^: 1+1 = -i+-i-. (25) 
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Introducing this result into eq. (l2^1) we get 



/= / dz{(l + (l-z) 2 )z^-(l + z 2 )(l-z)^ + \±^(l-zy-f q ( lq )5(l-z)}. (26) 

JO L — Z 

The integrals over the first two terms on the right-hand side (RHS) of (1261) exactly cancel 
as one sees by using the change of variable z — > 1 — z in one of them and the integral over 
the last two terms on the RHS of (|2"61) vanishes from the normalization in eq.([2]). Thus 
we conclude that 

/ = 0. (27) 
The quark momentum sum rule is indeed satisfied. 

Having improved the IR divergence properties of P qq (z) and PG q {z), we now turn to 
Pgg(z) and P q g(z). We first note that the standard formula for P q c(z), 

P q c(z) = \(z 2 + (1 - z)% (28) 

is already well-behaved (integrable) in the IR regime. Thus, we do not need to improve it 
here to make it integrable and we note that the singular contributions in the other kernels 
are expected to dominate the evolution effects in any case. We do not exclude improving 
it for the best precision [21] and we return to this point presently. 

This brings us then to Pgg{ z )- Its lowest order form is 

Pgg(z) = 2C G (— + -^- + z(1- z)) (29) 
z 1 — z 

which again exhibits unintegrable IR singularities at both z = 1 and z = 0. (Here, Cq 
is the gluon quadratic Casimir invariant, so that it is just N c = 3.) If we repeat the 
QCD exponentiation calculation carried-out above by using the color representation for 
the gluon rather than that for the quarks, i.e., if we apply the exponentiation analysis in 
the Appendix to the squared amplitude for the process G — > G(z) + G{1 — z), we get the 
exponentiated un-normalized result 

P GG (z) = 2C G F YFS {lG)e^ G (— Z ™ + -^-(1 - z)™ + V +7G (1 - z) + z(l - z) 1+ ^)) 

\ z 1 — z 2 / 

(30) 

wherein we obtain the jg and Sg from the expressions for 7 9 and 5 q by the substitution 
Cf — > Cq- 

1g = C G ^t = ^ (31) 

G = T^^ l y-^ (32) 

We see again that exponentiation has again made the singularities at z — 1 and z = 
integrable. 



C G 


—t = 

7T 


±C G 

A) 


7g 


a 

+ - 


sC G y 


2 


7T 1 3 
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To normalize P GG , we take into account the virtual corrections such that the gluon 
momentum sum rule 



/ dzz{2n f P qG {z)+P GG {z))=0 
Jo 



is satisfied. This gives us finally the IR-improved result 

P GG {z) = 2C G F YFS ( lG )e^ G {^^z^ + ^—(1 - zfe 

+ ^ 1+7G (1 -z) + z(l - ^) 1+7G ) - fahaW - z)} 

where for /g( 7 g) we get 

Tlf lr 2 1 



(33) 



(34) 



6C G F yF5 ( 7G ) 7g(1+7g)(2 + 7g) (1 + 7g)(2 + 7g) 

1 1 



(35) 



2(3 + 7 g)(4 + 7g) (2 + 7g)(3 + 7g)(4 + 7g) ' 

It is these improved results in eqs. (l2Q|22)34p for P qq (z), P Gq (z) and P GG (z) that we use 
together with the standard result in (129]) for P qG (z) as the IR-improved DGLAP theory. 

For clarity we summarize at this point the new IR-improved kernel set as follows: 

'l + z 2 



P%»{z) = C F F YFS { lq )e 



L 2 6 'l 



1 - z 



(1 - Z)T - f q { lq )5{\ - Z) 



(36) 



1 + (1 - zf 

z 

1 — z „,„ 



Pff(z) = C F F YFS { lq )e-^ ^ KX *> (37) 



Pgg(z) = 2C G F YFS { lG )e^{ z* + (1 



zV G 



+ \{z l+1G d -z) + z(l - ^) 1+7G ) - f G (lG)S(l - z)}, (38) 
P qG {z) = l -{z^{l-z)% (39) 

where we have introduced the superscript exp to denote the exponentiated results hence- 
forth. 

Returning now to the improvement of P qG (z), let us record it as well for the sake of 
completeness and of providing better precision. Applying eq.( l80|) to the process G —> q+q, 
we get the exponentiated result 

Pg(z) = F YFS { lG )e^ l -{z\l - z)*> + (1 - zfz^}. (40) 

The gluon momentum sum rule then gives the new normalization constant for the P G X G 
via the result 

7g(7g) = £r 



Cg(1 + 7g)(2 + 7g)(3 + 7g) 7g(1 + 7g)(2 + 7g) (1+7g)(2 + 7g) ^ 



1 1 



2(3 + 7 g) (4 + 7g) (2 + 7g) (3 + 7g) (4 + 7g) ' 
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The constant fa should be substituted for fa in Pqq whenever the exponentiated result 
in (HOl) is used. These results (1391) . (140j) . and (141j) are our new improved DGLAP kernel 
set, with the option exponentiating P q a as well. Let us now look into their effects on 
the moments of the structure functions by discussing the corresponding effects on the 
moments of the parton distributions. 

We know that moments of the kernels determine the exponents in the logarithmic 
variation [3,22] of the moments of the quark distributions and, thereby, of the moments 
of the structure functions themselves. To wit, in the non-singlet case, we have 

dM™(t) _ a,(t) A „ SKS{t) 



dt 2ir 
where 

M^ s (t)= ! dzz n - 1 q NS (z,t) (43) 



and the quantity A^ s is given by 

A% s = [ dzz n ~ l P e q ^{z) 
Jo 







= C F F YFS { lg )e^[B(n, lq ) + B{n + 2, lq ) - f q { lq )\ (44) 
where B(x,y) is the beta function given by 

B(x,y) = T(x)T(y)/T(x + y) 
. This should be compared to the un-IR-improved result [3,22]: 



1 1 A 1 

2 



2 n(n + l) 4-^ 



n 

3=2 J J 



j 



(45) 



The asymptotic behavior for large n is now very different, as the IR- improved expo- 
nent approaches a constant, a multiple of — f q , as we would expect as n — >• oo because 
lim^oo z n ~ 1 = for < z < 1 whereas, as it is well-known, the un-IR-improved result in 
(j45]) diverges as — 2CVln?7, as n — >• oo. The two results are also different at finite n: for 
n = 2 we get, for example, for a s = .118 [23], 

a ns _ \C F (-1.33) , un-IR-improved 
2 " |C F (-0.966) , IR-improved 

so that the effects we have calculated are important for all n in general. For completeness, 
we note that the solution to (I42p is given by the standard methods as 



t,t large with t»t \ Ot s {t) 
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where Ei(x) = dre r /r is the exponential integral function, 

jF -F YFS { lq )e^[B{n, lq ) + B{n + 2, 7ff ) - /,( 7 ,)] 



o 



with 



Ji (a B (t ) - a B (t)) 
2 ln(a,(t )/a s (*)) 

7T V 3 2 



(48) 



. We can compare with the un-IR- improved result in which the last line in eq. (|47I) holds 
exactly with a' n = 2A^ S ° /f] . Phenomenologically, for n = 2, taking Q = 2GeV and 
evolving to Q — lOOGeV, if we set Aq CD — -2GeV and use rif = 5 for definiteness of 
illustration, we see from eqs. (I47II48P that we get a shift of the respective evolved NS 
moment by ~ 5%, which is of some interest in view of the expected HERA precision [24]. 

We give now the remaining elements of the anomalous dimension matrix in its 'best' 
IR-improved form for completeness: 



n + 7g - 1 



B(3,7i + 7 ,-l) 



(49) 



An 9 = C dzz^P^z) = C F F YFS { lq )e^ 
Jo 

A G n G = [ dzz n - 1 P^(z)=2C G F YFS ( lG )e^{B(n + l, lG )+B(n + lG -l,2) 
Jo 

+ l -{B{n + 1, 7G + 2) + B(n + lG + 1, 2)) - f(<y G )}, (50) 

2n f Af = 2n f [ dzz n ~ l P™ p (z) = 2T{F)F YFS { lG )e^ G {B{n + 2, 1 + j G ) + B(n + 7g , 3)) 
Jo 

(51) 

where T(F) = |n/. We note that the un-exponentiated value of the last result in eq. floTT) 

is a well-known one [3,22], 2T{F) n (^Xi)(n+2) ■> an< ^ ^ would be used whenever we do not 
choose to exponentiate P qG . We will investigate the further implications of these IR- 
improved results for LHC physics elsewhere [21]. 

In the discussion so far, we have used the lowest order DGLAP kernel set to illustrate 
how important the resummation which we present here can be. In the literature [25,26], 
there are now exact results up to 0(a^) for the DGLAP kernels. The question naturally 
arises as to the relationship of our work to these fixed-order exact results. We stress first 
that we are presenting an improvement of the fixed-order results such that the singular 
pieces of the any exact fixed-order result, i.e., the parts, are exponentiated so that 

they are replaced with integrable functions proportional to (1 — z) 1 ^ 1 with 7 positive as 
we have illustrated above. Since the series of logs which we resum to accomplish this has 
the structure a e s t e \n n (l — z), £ > n these terms are not already present in the results 
in Refs. [25,26]. As we use the formula in eq. flHOj) . there will be no double counting if 
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we implement our IR-improvement of the exact fixed-order results in Refs. [25,26]. The 
detailed discussion of the application of our theory to the results in Refs. [25, 26] will 
appear elsewhere [21]. For reference, we note that the higher order kernel corrections in 
Refs. [25, 26] are perturbatively related to the leading order kernels, so one can expect 
that the size of the exponentiation effects illustrated above will only be perturbatively 
modified by the higher order kernel corrections, leaving the same qualitative behavior in 
general. 

In the interest of specificness, let us illustrate the IR-improvement of P qq when cal- 
culated to three loops using the results in Refs. [25,26]. Considering the non-singlet case 
for definiteness (a similar analysis holds for the singlet case) we write in the notation of 
the latter references 

oo 

P^s = PI + = EO n+lp - + ( 52 ) 



n=0 



where at order 0(a s ) we have 



P® + (z) = 2C F {-^±^- + - z)} (53) 

which shows that Pn° s )+ (z) agrees with the unexponentiated result in ([7j) for P qq except 
for an overall factor of 2. We use this latter identification to connect our work with that 
in Refs. [25,26] in the standard methodology. In Refs. [25,26], exact results are given for 
Pns + {z), and in Refs. [26] exact results are given for Pns + (z). When we apply the result 
in flBU]) to the squared amplitudes for the processes g — > q + X, q ^ q + X', we get the 
exponentiated result 

■^)2P^(z) + F YFS ( lq )e^ [(f^Ul " zy q Pi l s )+ (z) + B 2 6(l ~ *)} 
+ (^) 3 {(1 - zy«Pi? + (z) + B 3 5(l - z)}] 

where P^ v {z) is given in (J39l) and the resummed residuals Pn) + , i — 1,2 are related to 
the exact results for Pn) + , i = 1, 2, as follows: 

P$ + (z) = P^} + {z) - B 1+i S(l -z) + A% + (z) (55) 



./~\cm Ire T /■ ii" / 

ns , 



where 



A« + W = -AC^i 1 -^ - f q 5(l - z)} 



l + z 2 (56) 



and 



AE )+ (^) = -AC^fij^ - f q S(l - z)} 
-2tt5 1 P^ + (z) 

B 2 = B 2 + 4C F 7r5 1 f q 
B 3 = B 3 + ACp^ff, - 2tx8 x B 2 . 



(57) 
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Here, the constants B; t 



% = 2, 3 are given by the results in Refs. [25,26] as 



B 3 



B 2 




(58) 



where ( n is the Riemann zeta function evaluated at argument n. The detailed phe- 
nomenological consequences of the fully exponentiated 2- and 3-loop DGLAP kernel set 
will appear elsewhere [21]. 

In summary, we have used exact re- arrangement of the QCD Feynman series to isolate 
and resum the leading IR contributions to the physical processes that generate the evolu- 
tion kernels in DGLAP theory. In this way, we have obviated the need to employ artificial 
mathematical regularization of the attendant IR singularities as the theory's higher or- 
der corrections naturally tame these singularities. The resulting IR-improved anomalous 
dimension matrix behaves more physically for large n and receives significant effects at 
finite n from the exponentiation. 

We in principle can make contact with the moment-space resummation results in 
Ref. [27] but we stress that these results have necessarily been obtained after making a 
Mellin transform of the mathematical artifact which we address in this paper. Thus, the 
results in Ref. [27] do not in any way contradict the analysis in this paper. 

We note that the program of improvement of the hadron cross section calculations for 
LHC physics advanced herein should be distinguished from the results in Refs. [28,29]. 
Indeed, recalling the standard hadron cross section formula 



where {Fg(x)} are the respective parton densities and a(xiX2s) is the respective reduced 
hard parton cross section, the resummation results in Refs. [28,29] address, by summing 
the large logs in Mellin transform space, the x±x 2 —> 1 limit of <j{x\X2s) whereas the 
results above address the improvement, by resummation in x-space, of the calculation 
of the parton densities {Fi(x)} for all values of x. Thus, the program of improvement 




(59) 
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presented above is entirely complementary to that in Refs. [28, 29] and both programs of 
improvement are needed for precision LHC physics. 

Finally, we address the issue of the relationship between the re- arrangement that 
we have made of the exact leading-logs in the QCD perturbation theory and the usual 
treatment in the non-exponentiated DGLAP theory. If one expands out the exponentiated 
kernels, using the distribution identity 



one can see that for example P qq and P^ p agree to leading order, so that the leading log 
series which they generate for the respective NS structure functions also agree through 
leading order in — L where L is the respective big log in momentum-space. At higher 
orders then, we have a different result for the {Fi}, let us denote them by {F'i}, and a 
different result for the reduced cross section, let us denote it by a', such that we get the 
same perturbative QCD cross section, 



order by order in perturbation theory. The exponentiated kernels are used to factorize 
the mass singularities from the unfactorized reduced cross section and this generates 
a' instead of the usual a whose factorized form is generated using the usual DGLAP 
kernels. We thus have the same leading log series for a as does the usual calculation 
with un-exponentiated DGLAP kernels. We have an important advantage: the lack of 
+- functions in the generation of the configuration space functions {F'i, a'} means that 
these functions lend themselves more readily to Monte Carlo realization to arbitrarily soft 
radiative effects, both for the generation of the parton shower associated to the {F'i} and 
for the attendant remaining radiative effects in a'. 

Further consequences of our results for LHC physics will be presented elsewhere [21]. 
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Figure 2: The process Q'Q — > Q'" + Q" + n(G). The four-momenta are indicated in the 
standard manner: q 1 is the four-momentum of the incoming Q, q 2 is the four-momentum 
of the outgoing Q", etc., and Q — u,d, s, c, b, G. 

Appendix 



In this Appendix, we present the new QCD exponentiation theory which has been devel- 
oped in Refs. [16,17] as it is not generally familiar. The goal is to make the current paper 
self-contained. 

For definiteness, we will use the process in Fig. 2, Q'(pi)Q(qi) — > Q'" {P2)Q" {(It) + 
Gi(ki) ■ ■ ■ G n (k n ), as the proto-typical process, where we have written the kinematics 
as it is illustrated in the figure. This process, which dominates processes such as tt 
production at FNAL, contains all of the theoretical issues that we must face at the parton 
level to establish , as an extension of the original ideas of Yennie, Frautschi and Suura 
(YFS) [18], QCD soft exponentiation by MC methods - applicability to other related 
processes will be immediate. For reference, let us also note that, in what follows, we use 
the GPS conventions of Ref. [30] for spinors {u, v, u} and the attendant photon and gluon 
polarization vectors that follow therefrom: 



y/2 u- a (k)u a (P)' 



(62) 



with (3 2 = and £ defined in Ref. [30], so that all phase information is strictly known 
in our amplitudes. This means that, although we shall use the older EEX realization 
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of YFS MC exponentiation as defined in Ref. [31], the realization of our results via the 
the newer CEEX realization of YFS exponentiation in Ref. [31] is also possible and is in 
progress [21]. 

Specifically, the authors in Refs. [17] have analyzed how in the special case of Born 
level color exchange one applies the YFS theory to QCD by extending the respective YFS 
IR singularity analysis to QCD to all orders in a s . Here, unlike what was emphasized 
in Refs. [17], we focus on the YFS theory as a general re- arrangement of renormalized 
perturbation theory based on its IR behavior, just as the renormalization group is a general 
property of renormalized perturbation theory based on its UV(ultra-violet) behavior. We 
will thus keep our arguments entirely general from the outset, so that it will be immediate 
that our result applies to any renormalized perturbation theory in which the cross section 
under study is finite. 

Let the amplitude for the emission of n real gluons in our proto-typical subprocess, 
Q a + Q' a —> Q' n Q"'"' + n(G), where a, a, 7, and 7 are color indices, be represented by 

M^ a& = J2 M it & > ( 63 ) 

1 

is the contribution to Ai^ from Feynman diagrams with i virtual loops. Sym- 
metrization yields 

where this last equation defines symmetric function of its arguments arguments 

k\, ki. A will be our infrared gluon regulator mass for IR singularities; n-dimensional 
regularization of the 't Hooft-Veltman [32] type is also possible as we shall see. 

We now define the virtual IR emission factor SQCD{k) for a gluon of 4-momentum k, 
for the k — > regime of the respective 4-dimensional loop integration as in (164]) . such that 

limA; (p 7 ^x leading Casimir contribution 

S QCD (k)p% a& ) = 0, (65) 

where we have now introduced the restriction to the leading color Casimir terms at one- 
looplfl so that in the expression for the respective one-loop correction and in that for 
for SQCD{k) given in Refs. [17], only the terms proportional to Cp should be retained here 
as we focus on the ff— >ff case, where f denotes a fermion. (Henceforth, when we refer to 
k — ► gluons we are always referring for virtual gluons to the corresponding regime of 
the 4-dimensional loop integration in the computation of .) 

3 These correspond with maximally non-Abelian terms in Ref. [33] but computed exactly rather than 
in the eikonal approximation. 
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In Ref. [17], the respective authors have calculated <Sqcd(&) using the running quark 
masses to regulate its collinear mass singularities, for example; n- dimensional regulariza- 
tion of the 't Hooft-Veltman type is also possible for these mass singularities and we will 
also illustrate this presently. 

We stress that SQco{k) has a freedom in it corresponding to the fact that any function 
ASQCD{k) which has the property that lim^o k 2 ASQCD{k)p^ = may be added to it. 

Since the virtual gluons in pf^ are all on equal footing by the symmetry of this 
function, if we look at gluon £, for example, we may write , for k# — * (0, 0, 0, 0) = O while 
the remaining hi are fixed away from O, the representation 

pf 1 = S Q cD{h) * p^i(h, - • • , k t -i) + fi}(ki, ■■■ , kt-i; h) (66) 

where the residual amplitude P}(k\,--- , fo-i; ki) will now be taken as defined by this last 
equation. It has two nice properties: 

• it is symmetric in its first i — 1 arguments 

• the IR singularities for gluon i that are contained in Sqcd^i) are no longer con- 
tained in it. 

We do not at this point discuss the extent to which there are any further remaining 
IR singularities for gluon t in f5\{k\, ■ ■ ■ , kg_i, ki). In an Abelian gauge theory like QED, 
as has been shown by Yennie, Frautschi and Suura in Ref. [18], there would not be any 
further such singularities; for a non-Abelian gauge theory like QCD, this point requires 
further discussion and we will come back to this point presently. 

We rather now stress that if we apply the representation fl66|) again we may write 

= SQCD(ke)S QC D(ke-i) * pf} 2 {ki, ■ • ■ , h- 2 ) 
+SQCD(h)/3j_ l (k 1 , • • • , k e _ 2 ; k£-i) 
+SQCD(ki-i)f3}_ 1 (k 1 , ■■■ , ke_ 2 ; h) 

+0% (ki, • • • , k e _ 2 ; k(,_i, k e ), (67) 

where this last equation serves to define the function /3f • • • , k^_ 2 ; kg). It has two nice 
properties: 

• it is symmetric in its first I — 2 arguments and in its last two arguments kg—i, ki 

• the infrared singularities for gluons I — 1 and I that are contained in SQCD{ki-i) and 
SqcDiki) are no longer contained in it. 

Continuing in this way, with repeated application of (]66p . we get finally the rigorous, exact 
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rearrangement of the contributions to p\ as 

I 

pP = S QCD (ki) ■ ■ ■ S QC D(ki)P° + E II SQCD{kj)Pl(ki) + ■■■ 

i=i j^i 

+(3j(k 1 ,--- ,kt), (68) 
where the virtual gluon residuals Pl(k[, ■ ■ ■ , k^) have two nice properties: 

• they are symmetric functions of their arguments 

• they do not contain any of the IR singularities which are contained in the product 
Sqcd(K)-~ SQCD(ki). 

Henceforth, we denote /?• as the function for reasons of pedagogy. We can not stress 
too much that (|68p is an exact rearrangement of the contributions of the Feynman diagrams 

(n) 

which contribute to p) ; it involves no approximations. Here also we note that the question 
of the absolute convergence of these Feynman diagrams from the standpoint of constructive 
field theory remains open as usual. Yennie, Frautschi and Suura [18] have already stressed that 
Feynman diagrammatic perturbation theory is non-rigorous from this standpoint. What we do 
claim is that the relationship between the YFS expansion and the usual perturbative Feynman 
diagrammatic expansion is itself rigorous even though neither of the two expansions themselves 
is rigorous. 

Introducing (|68p into (163j) yields a representation similar to that of YFS, and we will call it 
a "YFS representation", 

oo 

M {n) =e a s B QCD J2m < f\ (69) 
3=0 

where we have defined 

/d 4 k 
m _ x 2 + ie) SQCD<Kk " > ^ 

and 

We say that ([69]) is similar to the respective result of Yennie, Frautschi and Suura in Ref. [18] and 
is not identical to it because we have not proved that the functions Pi(ki, k{) are completely 
free of virtual IR singularities. What have shown is that they do not contain the IR singularities 
in the product SQCD{k\) • • • SQCD{ki) so that does not contain the virtual IR divergences 
generated by this product when it is integrated over the respective 4j-dimensional j-virtual gluon 
phase space. In an Abelian gauge theory, there are no other possible virtual IR divergences; in 
the non-Abelian gauge theory that we treat here, such additional IR divergences are possible 
and are expected; but, the result (|69|) does have an improved IR divergence structure over (j63j) 
in that all of the IR singularities associated with SQCD^k) are explicitly removed from the sum 

over the virtual IR improved loop contributions to all orders in a s (Q). 
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Turning now to the analogous rearrangement of the real IR singularities in the differential 
cross section associated with the M^ n \ we first note that we may write this cross section as 
follows according to the standard methods 

e 2a s ReB QCD , » ^ « 

° = n\ ./ 11 ffc 2 , + A 2 ) 1 / 2 + qi ~ P2 ~ Q2 ~ ^ %) 



m=l v m ' " y i=l 

p k '{pi,qi,P2,q2,h, - ■ ■ ,kn) — fTTi — : (72) 



where we have defined 



oo 

p ( - n \ P l,qi,P2,q2,h,--- ,K)= Yl llE^H 2 ( 73 ) 

color, spin j=0 

in the incoming QQ' cms system and we have absorbed the remaining kinematical factors 
for the initial state flux, spin and color averages into the normalization of the amplitudes 
M {n) for reasons of pedagogy so that the p^ are averaged over initial spins and colors 
and summed over final spins and colors. We now proceed in complete analogy with the 
discussion of p™ above. 

Specifically, for the functions p^'ipi, ?i,£>2, <?2, ki, ■ ■ • , k n ) = p( n '(ki, ■ ■ ■ , k n ) which are 
symmetric functions of their arguments hi, ■ ■ ■ , k n , we define first, for n — 1, 

lim k 2 ( P^'(k) (leading Casimir contribution — SQCD{k)p^ J = 0, (74) 
|jfc|->0 V / 

where the real infrared function Sqod^) is rigorously defined by this last equation and is 
explicitly computed in Refs. [17], wherein we retain here only the terms proportional to 
Cf from the result in Ref. [17] ; like its virtual counterpart Sqct>(&) it has a freedom in 
it in that any function ASQCD{k) with the property that lim^^ k 2 AS QCD{k) = may 
be added to it without affecting the defining relation (174]) . 

We can again repeat the analogous arguments of Ref. [18], following the corresponding 
steps in fl66|) - fP7T]) above for Sqcd to get the "YFS-like" result 



CO 



_ SUM IR (QCD)Y" I f] ^h. f ^ P ^-(Pi+gl-P2-g2-Efcj)+gQGD 

hJ l\ "1 J < 2 ') 4 (75) 

/3n(^l? ■ • • j kn) n rT 
P2 <?2 



with 



SUM IR (QCD) = 2a s ReB QC D + 2a s B QCD (K max ) 
f d 3 k ~ 

2a s B QC D(K niax ) = / -^-SqcD(k)6(K maiX - k), 
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/rfik r i 

-£-Sqcd(AO [e~ iy - k - 0(K max - k)\ , (76) 

ifl, = ^(i-W) _ 2a s ReB QCD da B , 

I fa = da B1 - S QC iy(k)da B , ... (77) 

where the (5 n are the QCD hard gluon residuals defined above; they are the non-Abelian analogs 
of the hard photon residuals defined by YFS. Here, for illustration, we have recorded the relation- 
ship between the n , n = 0,1 through 0(a s ) and the exact one-loop and single bremsstrahlung 
cross sections, da^ 1 " 100 ^, da B1 , respectively, where the latter may be taken from Ref. [34] We 
stress two things about the right-hand side of ([75]) : 

• It does not depend on the dummy parameter K max which has been introduced for can- 
cellation of the infrared divergences in SUMjr(QCD) to all orders in a s (Q) where Q is 
the hard scale in the parton scattering process under study here. 

• Its analog can also be derived in our new CEEX [31] format. 

We now return to the property of (|75p that distinguishes it from the Abelian result derived 
by Yennie, Frautschi and Suura - namely, the fact that, owing to its non-Abelian gauge theory 
origins, it is in general expected that there are infrared divergences in the (5 n which were not 
removed into the Sqcd, Sqcd when these infrared functions were isolated in our derivation of 

sa. 



More precisely, the left-hand side of (1751) is the fundamental reduced parton cross 
section and it should be infrared finite or else the entire QCD parton model has to be 
abandoned. 

There is an observation in the literature [35] that unless we use the approximation 
of massless incoming quarks, the reduced parton cross section on the left-hand side of 
(1751) diverges in the infrared regime at 0(a 2 s (Q)). We do not go into this issue here but 
either use the quark masses strictly as collinear limit regulators so that they are set to 
zero in the numerators of all Feynman diagrams in such a way that the limit lim m 2/ £ ;2^ , 
where E q is the quark energy, is taken everywhere that it is finite or, alternatively, we 
use n-dimensional methods to regulate such divergences while setting the quark masses to 
zero as that is an excellent approximation for the light quarks at FNAL and LHC energies 
- we take this issue up elsewhere. 

From the infrared finiteness of the left-hand side of (175]) and the infrared fmiteness of 
SUMjr(QCD), it follows that the quantity 

Aft = p -sum ir (qcd) HfT 

must also be infrared finite to all orders in a s . 

As we assume the QCD theory makes sense in some neighborhood of the origin for 
a s , we conclude that each order in a s must make an infrared finite contribution to da cxp . 
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At 0{aP s {Q)) , the only contribution to da cxp is the respective Born cross section given 
by /9q in (1751) and it is obviously infrared finite, where we use henceforth the notation 
fin^ to denote the 0(a e s (Q)) part of (3 n . Thus, we conclude that the lowest hard gluon 
residual (3^ is infrared finite. 

Let us now define the left-over non-Abelian infrared divergence part of each contribu- 
tion Pn via 

m = $a> + m i) 

where the new function is now completely free of any infrared divergences and the 
function D/3 n ^ contains all left-over infrared divergences in (3$ which are of non-Abelian 
origin and is normalized to vanish in the Abelian limit / {, c — > where f a t, c are the group 
structure constants. 

Further, we define Dfln by a minimal subtraction of the respective IR divergences 
in it so that it only contains the actual pole and transcendental constants, 1/e — Ce for 
e = 2 — d/2, where d is the dimension of space-time, in dimensional regularization or In A 2 
in the gluon mass regularization. Here, Ce is Euler's constant. 

For definiteness, we write this out explicitly as follows: 

n+e 



dPh D^^Y.dfln'iX 2 



where the coefficient functions d™ ,e are independent of A for A — > and dPh is the 
respective n-gluon Lorentz invariant phase space. 

At 0(a^(Q)), the IR finiteness of the contribution to da exp then requires the contri- 
bution 

„ n 1 I „ ,3, n—i 1 l+i 

£=0 ' j = l Jk 3> K max J j = Q ' j=£+l 

d% s(n-t-i) ( , , ) d 3 p 2 d 3 q 2 



J 



P2 12 



to be finite. 

From this it follows that 



n-i -, e+i 



M58 = /E^n/ >lf t^-wEj n 

J l=Q ■ j = l J Kj>K max J i=0 • j—£ + l 



Z?C -^(W...,^)^^ (79) 



P 2 <? 2 



is finite. Since the integration region for the final particles is arbitrary, the independent 
powers of the IR regulator ln(A 2 ) in this last equation must give vanishing contributions. 
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— (g) 

This means that we can drop the Df3n from our result (J75J) because they do not make a 
net contribution to the final parton cross section a exp . We thus finally arrive at the new 
rigorous result 



s SUM m (QCD) 



oo „ n 

E/n 

71=0 ^ j 



<Pkj_ f d A y iy . {pi+qi _ P2 _ q2 _ TJe . )+Dam 



A kj J (2tt) 



(80) 



* 0n{kl, • • • j k n )— -t: -7j- 
P2 Q2 

where now the hard gluon residuals f3 n (ki, . . . , k n ) defined by 



Pn{ki , ■ ■ ■ , k n — ^ (3^ ^ (k\ , . . . , k n ) 



(81) 



are free of all infrared divergences to all orders in a s (Q). This is a basic result of this Appendix. 

We note here that, contrary to what was claimed in the Appendix of the first paper in 
Refs. [17] and consistent with what is explained in the third reference in [17], the arguments in 
the first paper in Refs. [17] are not sufficient to derive the respective analog of eq. (l80p ; for, they 
did not really expose the compensation between the left over genuine non-Abelian IR virtual and 
real singularities between J dPh/3 n and J dPhj3 n+ \ respectively that really distinguishes QCD 
from QED, where no such compensation occurs in the (5 n residuals for QED. 

We point-out that the general non-Abelian exponentiation of the eikonal cross sections in 
QCD has been proven formally in Ref. [33]. The contact between Ref. [33] and our result (|80p 
is that, in the language of Ref. [33], our exponential factor corresponds to the N=l term in the 
exponent of eq.(10) of the latter reference. One also sees immediately the fundamental difference 
between what we derive in (I80p and the eikonal formula in Ref. [33]: our result (I80p is an exact 
re-arrangement of the complete cross section whereas the result in eq.(10) of Ref. [33] is an 
approximation to the complete cross section in which all terms that could not be eikonalized 
and exponentiated have been dropped. 
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